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Abstract 

We propose a calculus of local equations over one-way computing pat- 
terns [RBB03j, which preserves interpretations, and allows the rewriting of 
any pattern to a standard form where entanglement is done first, then mea- 
surements, then local corrections. We infer from this that patterns with no 
dependencies, or using only Pauli measurements, can only realise unitaries 
belonging to the Clifford group. 



1 Introduction 

The one-way model centres on 1-qubit measurements as the main ingredi- 
ent of quantum computation |RBB03j . and is believed by physicists to lend 
itself to easier implementations |Nie04| IND041 IBR041 ICAJ04] . During com- 
putations, measurements and local corrections are allowed to depend on the 
outcomes of previous measurements. 

We first develop a notation for such classically correlated sequences of entan- 
glements, measurements, and local corrections. Computations are organised 
in patterns, and we give a careful treatment of pattern composition and ten- 
sor products (parallel composition) of patterns. We show next that such 

*This work was partially supported by the PREA, MITACS, ORDCF and CFI projects. 
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pattern combinations reflect the corresponding combinations of unitary op- 
erators. An easy proof of universality, based on a family of 2-qubit patterns 
follows. 

So far, this constitutes mostly a work of clarification of what was already 
known from the series of papers introducing and investigating the properties 
of the one-way model |RBB03j . However, we work here with an extended no- 
tion of pattern, where inputs and outputs may overlap in any way one wants 
them to, and this obtains more efficient - in the sense of fewer qubits - imple- 
mentations of unitaries. Specifically, our generating set consists of two simple 
patterns, each one using only 2 qubits. From it we obtain a 3 qubits realisa- 
tion of the R z rotations and a 14 qubit implementation for the controlled-?/ 
family: a very significant reduction over the known implementations. 
However, the main point of this paper is to introduce alongside our notation, 
a calculus of local equations over patterns that exploits the fact that 1 -qubit 
^-measurements are closed under conjugation by Pauli operators. We show 
that this calculus is sound in that it preserves the patterns interpretations. 
Most importantly, we derive from it a simple algorithm by which any general 
pattern can be put into a standard form where entanglement is done first, 
then measurements, then corrections. 

The consequences of the existence of such a procedure are far-reaching. First, 
since entangling comes first, one can prepare the entire entangled state needed 
during the computation right at the start: one never has to do "on the fly" 
entanglements. Second, since local corrections come last, only the output 
qubits will ever need corrections. Third, the rewriting of a pattern to stan- 
dard form reveals parallelism in the pattern computation. In a general pat- 
tern, one is forced to compute sequentially and obey strictly the command 
sequence, whereas after standardisation, the dependency structure is relaxed, 
resulting in low depth complexity. Last, the existence of a standard form for 
any pattern also has interesting corollaries beyond implementation and com- 
plexity matters, as it follows from it that patterns using no dependencies, 
or using only the restricted class of Pauli measurements, can only realise a 
unitary belonging to the Clifford group. 

Acknowledgements: Elham Kashefi wishes to express her gratitude to 
Quentin for letting her collaborate with his father, Vincent Danos, during 
their stay in Canada. Prakash Panangaden wishes to express his gratitude 
to EPSRC for supporting his stay in Oxford where this collaboration began. 
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2 Computation Patterns 



We first develop a notation for 1-qubit measurement based computations. 
The basic commands one can use are: 

— 1-qubit measurements M" 

— 2-qubit entanglement operators 

— and 1-qubit Pauli corrections X i , Z i 

The indices i, j represent the qubits on which each of these operations apply, 
and a is a parameter in [0, 2tt}. Sequences of such commands, together with 
two distinguished — possibly overlapping — sets of qubits corresponding to 
inputs and outputs, will be called measurement patterns, or simply patterns. 
These patterns can be combined by composition and tensor product. 
Importantly corrections and measurements are allowed to depend on previous 
measurement outcomes. We shall prove later that patterns without those 
classical dependencies can only realise unitaries that are in the Clifford group. 
Thus dependencies are crucial if one wants to define a universal computing 
model; that is to say a model where all finite-dimensional unitaries can be 
realised, and it is also crucial to develop a notation that will handle these 
dependencies gracefully. 



2.1 Commands 

The entanglement commands are defined as := AZy, while the correction 
commands are the Pauli operators X i and Z { . 

A 1-qubit measurement command, written M", is given by a pair of comple- 
ment orthogonal projections, on: 

l+«> := ^(|0) + e fa |l)) (1) 
K> := ^(|0)-e fa |l» (2) 

It is easily seen that \+ a ), \— a ) form an orthonormal basis in C 2 , so they 
indeed define a 1-qubit measurement (of rank 2 n_1 , if n is the number of 
qubits in the ambient computing space). Measurements here will always be 
understood as destructive measurements, that is to say the concerned qubit 
is consumed in the measurement operation. 

The outcome of a measurement done at qubit i will be denoted by Sj e Z 2 . 
Since one only deals with patterns where qubits are measured at most once 
(see condition (Dl) below), this is unambiguous. We take the convention 
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that Si — if under the corresponding measurement the state collapses to 
|+ Q ), and Si — 1 if to |— a ). 

Outcomes can be summed together resulting in expressions of the form s = 
Y^iiei Si wm ch we call signals, and where the summation is understood as 
being done is Z2. We define the domain of a signal as the set of qubits it 
depends on. 

Dependent corrections will be written X? and Z s { with s a signal. Their 
meaning is that X® = Zf = I (no correction is applied), while X\ = X i and 

Z} = z t . 

Dependent measurements will be written *[M"] S , where s and t are signals. 
Their meaning is as follows: 

*[Mf] s := Mt lYa+tn (3) 

As a result, before applying a measurement, one has to know first all the 
measurements outcomes occurring in the signals s, t, then one has to compute 
the parity of s and t, and maybe modify a to one of —a, a + ir and —a + n. 
One can easily compute that: 

XMtXi = Mr* (4) 
ZiM^Z, = M? +w (5) 

so that the actions correspond to conjugations of measurements under X 
and Z. We will refer to them as the X and Z-actions. Note that these two 
actions are commuting, since —a + 7r = —a — n up to 2tt, and hence the 
order in which one applies them doesn't matter. Should one use other local 
corrections, then one would have here instead the corresponding actions on 
measurement angles. As we will see later, relations (j3J) and (jSJ) are key to 
the propagation of dependent corrections, and to obtaining patterns in the 
standard entanglement, measurement, correction form. Since measurements 
considered here are destructive ones, the equations simplify to MfXi = Mr a y 
and MfZ t = M? +7T . 

Another point worth noticing is that the domain of the signals of a depen- 
dent command, be it a measurement or a correction, represents the set of 
measurements which one has to do before one can determine the actual value 
of the command. 

Finally we note that we could work with general 1-qubit measurements, in- 
stead of the class defined above, sometimes called :q/-measurements. All the 
developments would carry through nicely, but we have not found so far any 
compelling reason for this additional generality. 
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2.2 Patterns 



Definition 1 Patterns consists of three finite sets V , I, O, together with two 
infective maps i : I —> V and o : O — ► V and a finite sequence of commands 
A n . . . A\ applying to qubits in V . 

The set V is called the pattern computation space, and we write f)y for the 
associated quantum state space ®ie\/C 2 . To ease notation, we will forget 
altogether about the maps t and o, and write simply /, O instead of t(I) and 
o(0). Note however, that these maps are useful to define classical manipu- 
lations of the quantum states, such as permutations of the qubits. The sets 
/, O will be called respectively the pattern inputs and outputs, and we will 
write Sjj, and S)o for the associated quantum state spaces. The sequence 
A n . . . A\ will be called the pattern command sequence. 
To run a pattern, one prepares the input qubits in some input state ip G S)i, 
while the non-input qubits are all set in the |+) state, then the commands 
are executed in sequence, and finally the result of the pattern computation is 
some G Sjo- There might be qubits in the pattern, which are neither inputs 
nor outputs qubits, and are used as auxiliary qubits during the computation. 
Usually one tries to use as few of them as possible, since these participate to 
the space complexity of the computation. 

Note that one does not require inputs and outputs to be disjoint subsets of 
V . This seemingly innocuous additional flexibility is actually quite useful 
to give parsimonious implementations of unitaries |DKP04| . While the re- 
striction to disjoint inputs and outputs is unnecessary, it has been discussed 
whether more constrained patterns might be easier to realise physically. Re- 
cent work |HKB04[ IBR041 ICA,)04| however, seems to indicate they are not. 



Here is an example of a pattern implementing the Hadamard operator H: 

fi := ({1,2},{1},{2},X 2 S1 M 1 °E 12 ) 

What is this pattern doing ? The first qubit is prepared in some input 
state ip, and the second in state |+), then these are entangled to obtain 
AZi2(^i <S> |+)2)- Once this is done, the first qubit is measured in the |+), 
| — ) basis. Finally an X correction is applied on the output qubit, depending 
on the outcome of the measurement. We will do this calculation in detail 
later. 



2.3 Pattern combination 

We are interested now in how one can combine patterns into bigger ones. 
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The first way to combine patterns is by composing them. Two patterns ^ 
and m ay be composed if V\ fl V 2 = 0\ = I 2 . Note that provided that tyi 
has as many outputs as ^2 has inputs, by renaming the pattern qubits, one 
can always make them composable. 

Definition 2 The composite pattern ^2^Pi is defined as: 
-V :=V 1 UV 2 , I = h, O = 2 , 

— commands are concatenated. 

The other way of combining patterns is to tensor them. Two patterns %$i 
and ty$ 2 m ay be tensored if V\ fl V 2 = 0. Again one can always meet this 
condition by renaming qubits in a way that these sets are made disjoint. 

Definition 3 The tensor pattern ^2 <8> ^Pi is defined as: 

— V = V x U V 2 , I = h U I 2 , and O = O x U 2 , 

— commands are concatenated. 

Note that all unions above are disjoint. Note also that, in opposition to the 
composition case, commands from distinct patterns freely commute, since 
they apply to disjoint qubits and are independent of each other, so when we 
say that commands have to be concatenated, this is only for definiteness. 

2.4 Pattern conditions 

One might want to subject patterns to various conditions: 

(DO) no command depends an outcome not yet measured; 
(Dl) no command acts on a qubit already measured; 
(D2) a qubit i is measured if and only if i is not an output; 
(EMC) commands occur Es first, then Ms, then Cs. 

The reader might want to check that our example fj satisfies all of the above. 
It is routine to verify that these conditions are preserved under composition 
and tensor. Conditions (DO) and (Dl) ensure that a pattern can always be 
run meaningfully. Indeed if (DO) fails, then at some point of the computation, 
one will want to execute a command which depends on outcomes that are 
not known yet. Likewise, if (Dl) fails, one will try to apply a command 
on a qubit that has been consumed by a measurement (recall that we use 
destructive measurements). Condition (D2) is there to make sure that at the 
end of running the pattern, the state will belong to the output space S)o, 
i.e., that all non-output qubits, and only them, will have been consumed by 
a measurement when the computation ends. 

Starting now we will assume that all patterns satisfy the definiteness con- 
ditions (DO), (Dl) and (D2), and will designate by (D) the conjunction of 
these three conditions. 
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Condition (EMC) is of a completely different nature. Patterns not respecting 
it will be called wild. 

Later on, we will introduce the measurement calculus and show a simple 
rewriting procedure turning any given wild pattern into an equivalent one 
which is in (EMC) form. We call this procedure standardisation, and also 
say that a pattern meeting the (EMC) condition is standard. 
Before turning to this matter, we need a clean definition of what it means for 
a pattern to implement or to realise a unitary operator, together with a proof 
that the way one can combine patterns is reflected in their interpretations. 
This is key to our proof of universality. 



3 Computing a pattern 

Besides quantum states which are vectors in some Sjy, one needs a classical 
state recording the outcomes of the successive measurements one does in a 
pattern. So it is natural to define the computation state space as: 

S := [j VjW S) v xZY 

where V, W range over finite sets. In other words a computation state is a 
pair q, T, where q is a quantum state and T is a map from some W to the 
outcome space Z2. We call this classical component T an outcome map and 
denote by the unique map in Z 2 . 



3.1 Commands as actions 

We need a few notations. For any signal s and classical state T e Z^, such 
that the domain of s is included in W, we take sr to be the value of s given 
by the outcome map T. That is to say, if s = J2i s i> then s r := X^T(i) 
where the sum is taken in Z 2 . Also if T e Z^, and x G Z 2 , we define: 

Y[x/i\(i) = x, T[x/i](j) = T(j) for j ^ % 
which is a map in Z^ . 

We may now see each of our commands as acting on S. 







AZ i:j q, T 


g,r 






<z,r 






<z,r 


t [Af a] 8 


(+«rl 4 g,r[o/< 


<z,r 


'[M°] s 


(-arli?>r[l/< 



7 



where ar = (— l) Sr a + trii following equation Q, and (^Ij is the linear form 
associated to ip applied at qubit i. Suppose q G fiy, for the above relations 
to be defined, one needs the indices i, j on which the various command apply 
to be in V. One also needs T to contain the domains of s and t, so that sr 
and tr are well-defined. This will always be the case during the run of a 
pattern because of condition (D). 

All commands except measurements are deterministic and only modify the 
quantum part of the state. The measurements actions on S are not deter- 
ministic, so that these are actually binary relations on S, and modify both 
the quantum and classical parts of the state. The usual convention has it 
that when one does a measurement the resulting state is renormalised, but 
we don't adhere to it here, the reason being that this way, the probability of 
reaching a given state can be read off its norm, and the overall treatment is 
simpler. 

We introduce an additional command called shifting: 

q,T g,r[r(f)+5r/i] 

It consists in shifting the measurement outcome at % by the amount sp. Note 
that the Z-action leaves measurements globally invariant, in the sense that 
\+a+ir), \— a+n) = \~a), \+a) ■ Thus changing a to a + 7i amounts to swap the 
outcomes of the measurements, and one has: 

'[MtY = St[M?] s (6) 

and shifting allows to split the t action of a measurement, resulting sometimes 
in convenient optimisations of standard forms. 



3.2 Computation branches 

Let ^3 be a pattern with computation space V, inputs /, outputs O and 
command sequence A n . . . A%. A complete pattern computation starts with 
some input state q in fij, together with the empty outcome map 0. The 
input state q is then tensored with as many |+)s as there are non-inputs in 
V, so as to obtain a state in the full space Sjy- Then commands in *P are 
applied in sequence. We can summarise the situation as follows: 



53j 



o 



prep 



v\o 



S 



To make this precise, say there is a ^-branch from g £ Sjj to g' £ i^o, written 
g ^<p g', if there is a sequence (g^, Tj) with 1 < i < n + 1, such that: 

q(g) |+...+),0 = g 1 ,r 1 

and for all i < n : gj, Tj — ^ g i+ i, 

thus — ><p is a binary relation on Sji x ^o- That it is a relation and not a 
map reflects the fact that measurements a priori introduce non determinism 
in the evolution of the quantum states. 

Specifically, if k is the number of measurements in ^3 (or equivalently the 
number of non-outputs qubits), there are at most 2 k branches in any given 
computation, and therefore a given q £ Sji is in relation with at most 2 k 
distinct q' £ S)o- The probability of a branch is defined to be ||g'|| 2 /lkll 2 (g 
being always assumed to be non zero). Indeed one has: 

£ UsT = imi 2 (7) 

{q'\q^yq'} 

since any action is either a unitary, thus a norm-preserving action, or a 
measurement which introduces a branching, and then if q projects to go and 
gi, under some M", ||g|| 2 = ||go|| 2 + ||gi|| 2 , so that the relation above is always 
preserved. 

Definition 4 One says the pattern is deterministic if for all q £ S^i, q' 
and q" £ S^o, whenever q -^*p g' and q — q" , then q' and q" only differ up 
to a scalar. 

Note that even when ^3 is deterministic, all branches might not be equally 
likely. When ty$ is deterministic, one defines a norm-preserving map £Ap from 
f)/ to S) by: 

U v (q) := ( 8 ) 

Note that when q — ><p g', g' 7^ 0, so that the definition above always make 
sense. Note also that because is deterministic, this map depends on the 
choice of q' only up to a global phase. One can further comment that since 
we took the convention not to renormalise measurement results, we have to 
do here a global renormalisation to define the pattern interpretation. 
One says that a deterministic pattern *P realises or implements Uy, or equiv- 
alently that £/fp is the interpretation of ^3. 
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This map U<g must actually be a unitary embedding, since all quantum de- 
finable deterministic transformations are unitaries. If a precise argument is 
needed here, one can rephrase all the definitions given so far in the language 
of density operators and completely-positive maps (cp-maps). Then a deter- 
ministic pattern will implement a cp-map preserving pure density operators. 
From the Kraus representation theorem for cp-maps, it is easy to see that 
such cp-maps are liftings of unitary embeddings. 



3.3 Short examples 

First we give a quick example of a deterministic pattern that has branches 
with different probabilities. The state space is {1,2}, with 1 = = {1}, 
while the command sequence is M^. Therefore, starting with input q, one 
gets two branches: 

Ma { !(l + e--)g,0[O/2] 
g®|+),0 -A I 

{ |(l-e--)g,0[l/2] 

Thus this pattern is indeed deterministic, and implements the identity up 
to a global phase, and yet the two branches have respective probabilities 
(1 + cosa)/2 and (1 — cosa)/2, which are not equal in general. 
Next, we return to the pattern Sj which we already took as an example. Let 
us consider for a start the pattern with same space {1,2}, same inputs and 
outputs / = {1}, O = {2}, and shorter command sequence M®E 12 . Starting 
with input q = (a|0)+6|l))|+), one has two computation branches, branching 
at M°: 

(a|O) + 6|l))|+),0 ^(a|OO) + a|Ol)+6|lO)-6|ll)),0 

M o j |((a + 6)|O> + (a-6)|l»,0[O/O] 
\ i((a-6)|O> + (a + 6)|l»,0[l/O] 

and since \\a + b\\ 2 + \\a — b\\ 2 = 2(||a|| 2 + \\b\\ 2 ), both transitions happen 
with equal probabilities \. Both branches end up with different outputs, so 
the pattern is not deterministic. However, if one applies the local correction 
X 2 on either of the branches ends, both outputs will be made to coincide. 
Let us choose to let the correction bear on the second branch, obtaining the 
example ft which we defined already. We have just proved H = 11$, that is 
to say 9) realises the Hadamard operator. 
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With our definitions in place, we first infer that pattern combinations corre- 
spond to combinations of their interpretations. From this an easy structured 
argument - that uses surprisingly simple patterns - for universality will follow. 

3.4 Composing, Tensoring and Interpretation 

Recall that two patterns ^3i, ^2 may be combined by composition provided 
tyi have as many outputs as has inputs. Suppose this is the case, and 
suppose further that and ^2 respectively realise some unitaries U\ and 
U 2 , then the composite pattern ^2^Pi realises U^U-y. 
Indeed, the two diagrams representing branches in tyi and ^2: 

fih - #Oi - #o 2 

can be pasted together, since 0\ = I 2 , and $)o 1 = fti 2 - But then, it is 
enough to notice 1) that preparation steps p 2 in commute to all actions 
in since they apply on disjoint sets of qubits, and 2) that no action taken 
in depends on the measurements outcomes in tyy. It follows that the 
pasted diagram describes the same branches as does the one associated to 
the composite ^2^1 • 

A similar argument applies to the case of a tensor combination, and one has 
that «p 2 ® ^Pi realises U 2 ®Ui. 

The same holds even for non-deterministic patterns considered as implement- 
ing cp-maps. But we will not be concerned with this generalised setting in 
this paper. 

4 Universality 

Consider the two following patterns: 

Z{a) := X^M^ a E 12 (9) 
A3 := E 12 (10) 

In the first pattern 1 is the only input and 2 is the only output, while in the 
second both 1 and 2 are inputs and outputs. Note that here we are taking 
advantage of allowing patterns with overlapping inputs and outputs. 

Proposition 5 The patterns 3(a) and A3 are universal. 
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First, we claim 3(a) and A3 respectively realise J (a) and AZ, with: 



J (a) 




) 



We have already seen in our example that 3(0) = Sj implements H = J(0), 
thus we already know this in the particular case where a = 0. The general 
case follows by the same kind of computation. The case of AZ is obvious. 
Second, we know that these unitaries form a universal set |DKP04j . There- 
fore, from the preceding section, we infer that combining the corresponding 
patterns will generate patterns realising all finite-dimensional unitaries. □ 
These patterns are indeed among the simplest possible. As a consequence, in 
the section devoted to examples, we will find that our implementations have 
often little space complexity. 

Remarkably, in our set of generators, one finds a single dependency, which 
occurs in the correction phase of 3(a). No set of patterns without any mea- 
surement could be a generating set, since such patterns can only implement 
unitaries in the Clifford group. Dependencies are also needed for universality, 
but we have to wait for the development of the measurement calculus in the 
next section to give a proof of this fact. 

5 The measurement calculus 

We turn to the next important matter of the paper, namely standardisation. 
The idea is quite simple. It is enough to provide local pattern rewriting rules 
pushing Es to the beginning of the pattern, and Cs to the end. 

5.1 The equations 

A first set of equations give means to propagate local Pauli corrections 
through the entangling operator E^. Because Ey = Eji, there are only 
two cases to consider: 



These equations are easy to verify and are natural since E^ belongs to the 
Clifford group, and therefore maps under conjugation the Pauli group to 
itself. 






(12) 



12 



A second set of equations give means to push corrections through measure- 
ments acting on the same qubit. Again there are two cases: 

*[Mf] s A7 = *[Mf] s+r (13) 
*[Mf] s Z[ = t+r [M?] s (14) 

These equations follow easily from equations (J3J) and (J3J). They express the 
fact that the measurements M" are closed under conjugation by the Pauli 
group, very much like equations (fTTf and (fT2j) express the fact that the Pauli 
group is closed under conjugation by the entanglements Eij. 
Define the following convenient abbreviations: 

[Mf] s := °[Mf] s , *[Mf] := '[Mf] , Mf := °[Mf]°, 
Mf := MP, Mf := M} 

Particular cases of the equations above are: 

MfX? = M* 

MfXf = [MfY = *[Mf] = M\Zl 

The first equation, follows from —0 = 0, so the X action on Mf is trivial; the 
middle equation, second row, is because — f is equal | + ix modulo 2ti, and 
therefore the X and Z actions coincide on Mf . So we obtain the following: 

*[Mf] s = '[Mf] (15) 
*[Mf] s = s+t [Mf] (16) 

which we will use later to prove that patterns with measurements of the form 
M x and M y may only realise unitaries in the Clifford group. 



5.2 The rewriting rules 

We now define a set of rewrite rules, obtained by directing the equations 
above: 

EijXf =>- X-ZjEij EX 

EijZf =>- Z-Eij EZ 

'[Mf] s Xr '[Mf] s+r MX 

^{MfYZl r+t [M«Y MZ 

to which we need to add the free commutation rules, obtained when com- 
mands operate on disjoint sets of qubits: 

EijA % => with A ^ E 

A%X? => XfA % with A ^ C 
A % Zf Z?A % with A 
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where k represent the qubits acted upon by command A, and are supposed 
to be distinct from i and j. 

Condition (D) is easily seen to be preserved under rewriting. 
Under rewriting, the computation space, inputs and outputs remain the 
same, and so are the entanglement commands. Measurements might be mod- 
ified, but there is still the same number of them, and they are still acting on 
the same qubits. The only induced modifications concern local corrections 
and dependencies. We also take due note that none of these equations may 
create dependencies. 



5.3 Standardisation 

Write ^3 =>- ^3', respectively ^3 ^3', if both patterns have the same type, 
and one obtains ^3"s command sequence from *}3's one by applying one, 
respectively any number, of the rules above. Say ^3 is standard if for no ^3', 

«p => qy. 

Because all our equations are sound, one has that whenever *}3 =^>* ^3', and 
both patterns are deterministic, then {Tip = Uy. 

One can show by a standard rewriting theory argument, that for all *}3, there 
exists a unique standard ^3', such that *}3 =^>* ^3', and moreover ^3' satisfies 
the (EMC) condition. Reaching the standard form takes at most quadratic 
time in the number of instructions in ^3. Details are given in the appendix. 



5.4 Signal shifting 

One can extend the calculus to include the shifting command Sf. This allows 
one to dispose of dependencies induced by the Z-action, and obtain some- 
times standard patterns with smaller depth complexity, as we will see in the 
next section devoted to examples. 

*m- =► si[M:y (17) 

X^S* slxf t+s * ls * ] (18) 
Z]S] > Slzf t+Si/Si] (19) 
*[Mj*] s S7 => si t[r+Si/si] [M^] slr+Si/si] (20) 

where s[t/si] is the substitution of Sj with t in s, s, t being signals. The 
first additional rewrite rule was already introduced as equation (JSJ), while 
the other ones are merely propagating the signal shift. Clearly also, one can 
dispose of Sj when it hits the end of the pattern command sequence. We will 
refer to this new set of rules as =^#. 
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6 Examples 



In this section we develop some examples illustrating both pattern compo- 
sition, pattern standardisation, and signal shifting. We compare our imple- 
mentations with the implementations given in the reference paper |RBB03j . 
To combine patterns one needs to rename their qubits as we already noticed. 
We use the following concrete notation: if *p is a pattern over {1, . . . , n}, and 
/ is an injection, we write . . . , f{n)) for the same pattern with qubits 

renamed according to /. We also write ?fi 2 ° ^Pi for pattern composition to 
ease reading. 



Teleportation. 

Consider the composite pattern 3(/3)(2, 3)o-ji(a)(l, 2) with computation space 
{1, 2, 3}, inputs {1}, and outputs {3}. We run our standardisation procedure 
so as to obtain an equivalent standard pattern: 

3(/3)(2,3)o3(a)(l,2) = X^M- (3 E 23 X^M{ a E 12 

=^>ex X 3 S2 M 2 ^X 2 1 Z^ 1 M 1 a E 23 E\ 2 
=> MX X^Z^[M^M^ a E 23 E 12 

Let us call the pattern just obtained 3(a, (3). If we take as a special case 
a = (3 = 0, we get: 

X 3 S2 Z S 3 1 M$M?E 23 E 12 

and since we know that 3(0) implements H and H 2 = I, we conclude that 
this pattern implements the identity, or in other words it teleports qubit 1 
to qubit 3. As it happens, this pattern obtained by self-composition, is the 
same as the one given in the reference paper |RBB03[ p. 14]. 



x-rotation. 

Here is the reference implementation of an x-rotation |RBB03| p. 17], R x (a): 



X^Z S ^ [M 2 - a } s -M*E 23 E 12 (21) 

with computation space V = {1, 2, 3}, {1}, {3}. There is a natural question 
which me might call the recognition problem, namely how do we know this is 
implementing R x (a) ? Of course there is the brute force answer to that, which 
we applied to compute our simpler patterns, and which consists in computing 
down all the four possible branches generated by the measurements at 1 
and 2. Another possibility is to use the stabiliser formalism as explained 
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in the reference paper |RBB03j . Yet another possibility is to use pattern 
composition, as we did before, and this is what we are going to do. 
We know that R x {at) = J(a)H up to a global phase, hence the composite 
pattern 3(a)(2, 3) o 55(1, 2) implements R x (a). Now we may standardise it: 

3(a)(2,3)o£(l,2) = X?M; a E 23 X?M?E 12 

^ EX XI 2 Zl 1 M^ a Xl 1 M^E 23 E 12 
^mx X?Z s ^[M^MfE 23 E 12 

obtaining exactly the implementation we started with. Since our calculus is 
preserving interpretations, we deduce that the implementation is correct. 

^-rotation. 

Now, we have a method here for synthesising further implementations, which 
we can use fir instance with another rotation R z (a). Again we know that 
R z (cx) = HR x (a)H, and we already know how to implement both compo- 
nents H and R x (a). 

Starting with the pattern 55(4, 5) o 9^ x (o;)(2, 3, 4) o 55(1, 2) we get: 

55(4,5)o91 c ( a )(2,3,4)o55(l,2) = 

55(4, 5)X?Z? [M«] l + S *MZE 34 E 23 X^M*E 12 ^ EX 

55(4, 5)X?Z? [M^M-X^E U Z S ^M^E 123 ^ EZ 

55(4, 5)Xf Z? [M«} 1+s *Z s 3 1 M 2 rX 2 sl M?E 12U ^ MX 

55(4, 5)Xf [M 3 a ] 1+S2 Z 3 Sl Mf Mf £ 1234 ^mz 

X^MtE^XfZ 8 ^ [M 3 <f +**Mf Mf E 12U ^ EX 

Xl*Zl*MZXl*Z?* [M«} 1+S *M*M?E 12U5 ^ MX 

X?Z? [MZY*Z?* [M«} l + S *M%M?E 12U5 ^ MZ 

X S4 Z s 3S2 [ M x]53*i [M 3 a ] 1+S2 Mf Mf E12345 

To ease reading E 23 E 12 is shortened to -E123, E 12 E 23 E U to £'1234, and *[Mf] 1+s 
is used as shorthand for t [Mf a ] s . 

Here for the first time, we see MZ rewritings, inducing the Z-action on 
measurements. The obtained standardised pattern can therefore be rewritten 
further using the extended calculus: 

X° 4 Z° 3S2 [M%] s ^ [M«] 1+S2 Mf Mf £ 1234 5 =^s 
X S 5 2+S4 Z° 1+S3 M£ [M^] 1+S2 MfMfEi2345 

obtaining again the pattern given in the reference paper |RBB03l p. 5]. 
However, just as in the case of the R x rotation, we also have R z (a) = HJ(a) 
up to a global phase, hence the pattern 55(2, 3)3(a)(l, 2) also implements 
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R z (a), and we may standardize it: 



55(2,3) o3(a)(l, 2) = X s 3 2 M^E 23 X^M{ a E 12 
^ EX X^Z^M^X^M^E^ 
^ MX X s 3 *Z^M*M^ a E 123 



obtaining a 3 qubits standard pattern for the ^-rotation, which is simpler 
than the preceding one, because it is based on the 3(a) generators. Since the 
z-rotation R z (a) is the same as the phase operator: 



up to a global phase, we also obtain with the same pattern an implementation 
of the phase operator. In particular, if a — |, using the extended calculus, 
we get the following pattern for P(|): X s 3 2 Z s 3 1+1 M% M\ E 123 . 

General rotation. 

The realisation of a general rotation based on the Euler decomposition of 
rotations as -R x (7)i? 2 (/3)i? x (a), would results in a 7 qubits pattern. We get 
a 5 qubits implementation based on the J (a) decomposition |DKP04| : 



3(0) (4, 5)3(a)(3, 4)303) (2, 3)3(t)(1, 2) 

X^MlE A ,X^M^E u X s ^M^E 23 X^MlE l2 ^ EX 

X^M^E 45 X S 4 3 M^E U X^M^X^Z^M^E 123 ^ MX 

X S 5 *M°E 45 X S /M°E U X S 3 *Z^ [M^Y 1 MJ Ei2 3 ^ EXZ 

XfMZEtsXFMZXpZfZ? [M^M^E 12U ^ MXZ 

X^MlE^Z^ [M?] s * [M^M^E 1234 ^ EXZ 

X^M2X s 4 3 Zl 2 Z^[M^ [M^M^E l2U5 ^ MXZ 

XI* Z^" [M 4 °] S1 [Mg] s > [M^M^E 12U5 ^ s 
X 5 S2+S4 Z| 1+S3 M° [M 3 Q ] S2 [M^MjE 1234b 



CNOT (AX). 

This is our first example with two inputs and two outputs. We use here the 
trivial pattern 3 with computation space {1}, inputs {1}, outputs {1}, and 
empty command sequence, which implements the identity over Sj 4 . 




R(a,f3,>y) 



J(0)J(a)J(P)J(j) 



The extended standardization procedure yields: 
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One has AX = (I <8> H) t\Z(l ® H), so we get a pattern using 4 qubits over 
{1,2,3,4}, with inputs {1,2}, and outputs {1,4}, where one notices that 
inputs and outputs intersect on the control qubit {1}: 

(2f(l)®fj(3,4))A3(l,3)p(l)®f)(2,3)) = XfMf E U E 13 X S 3 *M$E 23 

By standardising: 

X S 4 3 M£E U E 13 X S 3 >M*E 23 ^ EX 
X S /Z S ^M^E U X^M^E 13 E 23 ^ EX 
X^ Z^ Z s ^ M^X^ M 2 E 13 E 23 E 34 =^mx 
X s / Z? Z? M*M*E 13 E 23 E U 

Note that we are not using here the Ei 234 abbreviation, because the un- 
derlying structure of entanglement is not a chain. This pattern was already 
described in Aliferis and Leung's paper |AL04| . In their original presentation 
the authors actually use an explicit identity pattern (using the teleportation 
pattern 3(0, 0) presented above), but we know from the careful presentation 
of composition that this is not necessary. 



GHZ. 

We present now a family of patterns preparing the GHZ entangled states 
|0. . .0) + |1 . . . 1). One has: 

GHZ(n) = (H n AZ„_i„ ...H 2 f\Z 12 )\ + . . .+) 

and by combining the patterns for AZ and H, we obtain a pattern with 
computation space {1, 2,2', ... , n, n'}, no inputs, outputs {1, 2' , . . . , n'}, and 
the following command sequence: 

X s n ?M*E nn ,E 

(n— l)'n • • • X^M 2 E 22 iE\ 2 

Under that form, the only apparent way to run the pattern is to execute all 
commands in sequence. The situation changes completely, when we bring 
the pattern to extended standard form: 

XrfM%E nn rE( n _iy n . . . X 3 fM 3 E 33 iE 2 i 3 X 2 s fM 2 E 22 iEi 2 =>■ 
X^X%M*E nn iE{ n _\y n . . . X^ M 3 Z3 2 Mf E 3V E 2 i 3 E 22 > E 12 =^ 
X s jX s 2 ?MZE nn ,E {n _ iyn . . . X s 3 f S2 [MZ\M%E 33 ,E 2 , 3 E 22 ,E 12 ^ 
X S J... X 3 s , 3 X 2 r«-i [M*] ... 52 [Mf ] M*E nn ,E {nr . iyn . . . E 33 ,E 2I3 E 22 ,E X2 ^ s 
X^f +S3+ +Sn . . . Xy +S3 X2?M* . . . M 3 M 2 E nn >E{n-iy n . . . E 33 iE 2 i 3 E 22 iE 12 

All measurements are now independent of each other, it is therefore possible 
after the entanglement phase, to do all of them in one round, and in a 
subsequent round to do all local corrections. In other words, the obtained 
pattern has constant depth complexity 2. 
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Controlled-^/ 



This final example presents another instance where standardization obtains 
a low depth complexity. For any 1-qubit unitary U, one has the following 
decomposition of AU in terms of the generators J(a) |DKP04| : 

AC/12 = ^1^1 J2J2 '^2 2 '^2 2 '^2 AZ12J2 "2 ^2 2 '^2 A-Zl2<-^2 2 

with a' = a + /3+ ^ +l5 . By translating each J operator to its corresponding 
pattern, we get the following wild pattern for AU : 

X s c B M°zE BC X s B A Mx a ' E AB X s k > MfEjkX^Mr^Eij 
Xt h M| E hi X s h 9 MJ E gh Xg f M° f E fg E A f X s f M~ 1 £ e/ 

_ 7 7T + l5 + /3 /3 — ^ + 7T 

X^MP E de X s /M~ E cd X^M^E bc E Ab X-^M~ E ab 

Figure ^ shows the underlying entanglement graph for the AU pattern. In 
order to run the wild form of the pattern one needs to follow the graph struc- 
ture and hence one has to perform the measurement commands in sequence. 
Extended standardisation yields: 




O k 



Figure 1: The underlying entanglement graph for the AU pattern. 



rySi + Sg+Se + Sc + Sa V Sj+S h +S f +S d + S b V S B /ySA + Se+Sc 

M B M^ a 'M° [M^~ 7T ] Sh+s f +Sd+Sb ] i M^} S3+Se+Sc+Sa [M]} s s +Sd+Sb 

77 7 7T — 5 — {3 — P -j- <5 -j- 77 

M° f [M e "] Si+St [M]] s - +S - [M c 2 } Sb M^M a 2 
EBcEABEjkEijEhiE g hEfgEAfE e fEdeE c dE bc E ab EAb 

Figure El shows the dependency structure of the resulting standard pattern 
for AU, and one sees it has depth complexity 7. 

7 The no dependency theorems 

From standardization we can also infer results related dependencies. We start 
with a simple observation which is a direct consequence of standardisation. 
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Figure 2: The dependency graph for the standard AU pattern. 



Lemma 6 Let ^ be a pattern implementing some unitary U, and suppose 
's command sequence has measurements only of the M x and M y kind, then 
U has a standard implementation, having only independent measurements, 
all being of the M x and M y kind (therefore of depth complexity at most 2). 

Write ^P' for the standard pattern associated to ^3. By equations (jl5j) and 
(jl6|) . the X-actions can be eliminated from ^3', and then Z-actions can be 
eliminated by using the extended calculus. The final pattern still implements 
U, has no longer any dependent measurements, and has therefore depth 
complexity at most 2. □ 

Theorem 1 Let U be a unitary operator, then U is in the Clifford group iff 
there exists a pattern *}3 implementing U , having measurements only of the 
M x and M y kind. 

The "only if" direction is easy, since we have seen in the example section, 
standard patterns for AX, H and P (|) which had only M x and M y measure- 
ments. Hence any Clifford operator can be implemented by a combination of 
these patterns. By the lemma above, we know we can actually choose these 
patterns to be standard. 

For the "if" direction, we prove that U belongs to the normaliser of the Pauli 
group, and hence by definition to the Clifford group. In order to do so we use 
the standard form of ^3 written as ^3' = C^M^E^/ which still implements 
U, and has only M x and M y measurements. 

Let i be an input qubit, and consider the pattern ^3" = ^3Cj, where C, is 
either Xj or Zj. Clearly *}3" implements UCi. First, one has: 

for some non-dependent sequence of corrections C, which, up to free commu- 
tations can be written uniquely as C' C" , where C applies on output qubits, 
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and therefore commutes to M<p/ , and C" applies on non-output qubits (which 
are therefore all measured in M«p/). So, by commuting C both through M«p/ 
and C«p/ (up to a global phase), one gets: 

Cqy MswC E?fii C'qC?qi MsqiC" Esq/ 

Using equations (fTKj) . (fT6|) . and the extended calculus to eliminate the re- 
maining Z-actions, one gets: 

M V C ^ MCiS SM V , 

for some product S = YlueJ} °^ cons tant shiftings, applying to some 
subset J of the non-output qubits. So: 

CqCsw Mm'C E<w CqCvq' S M<gr Esgi 

where Cq is a further constant correction obtained by shifting Csp/ with S. 
This proves that ^3" also implements C'qCqU, and therefore UCi = C'qC'qU 
which completes the proof, since C' Cq is a non dependent correction. □ 
The only if part of this theorem already appears in previous work |RBB03t 
p.18]. 

We can further prove that dependencies are crucial for the universality of 
the model. Observe first that if a pattern has no measurements, and hence 
no dependencies, then it follows from (D2) that V — O, i.e., all qubits are 
outputs. Therefore computation steps involve only X, Z and AZ, and it is 
not surprising that they compute a unitary which is in the Clifford group. 
The general argument essentially consists in showing that when there are 
measurements, but still no dependencies, then the measurements are playing 
no part in the result. 

Theorem 2 Let ^ be a pattern implementing some unitary U, and sup- 
pose ^3 's command sequence doesn't have any dependencies, then U is in the 
Clifford group. 

Write ^3' for the standard pattern associated to ^3. Since rewriting is sound, 
still implements U, and since rewriting never creates any dependency, it 
still has no dependencies. In particular, the corrections one finds at the end 
of call them C, bear no dependencies. Erasing them off ^3', results in 
a pattern which is still standard, still deterministic, and implementing 
U' := C*U. 

Now how does the pattern ty" run on some input <f ? First <fi <g> | + . . . +} 
goes by the entanglement phase to some ip G fjy, and is then subjected to 
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a sequence of independent 1-qubit measurements. Pick a basis B spanning 
the Hilbert space generated by the non-output qubits $Jv\o an d associated 
to this sequence of measurements. 

Since $j v = #0 ® fiv\o and fiv\o = ©^esf^fe], where [4> b ] is the linear 
subspace generated by </>&, by distributivity, if> uniquely decomposes as: 

i> = E^eB & ® x & 

where 06 ranges over B, and e i}o- Now since ^3" is deterministic, there 
exists an x, and scalars A& such that a;& = A&a;. Therefore if> can be written 
if}'®x, for some ■0'. It follows in particular that the output of the computation 
will still be x (up to a scalar), no matter what the actual measurements are. 
One can therefore choose them to be all of the M x kind, and by the preceding 
theorem U' is in the Clifford group, and so is U — CU', since C is a Pauli 
operator. □ 

From this section, we conclude in particular that any universal set of patterns 
has to include dependencies (by the preceding theorem), and also needs to 
use measurements M a where a^O modulo | (by the theorem before). This 
is indeed the case for the universal set 3(a) and A3- 

8 Conclusion 

We presented a calculus for 1-qubit measurement based quantum comput- 
ing. We have seen that pattern combinations allow for a structured proof of 
universality, which also results in parsimonious implementations. We have 
shown further that our calculus defines a quadratic-time standardisation al- 
gorithm transforming any pattern to a standard form where entanglement 
is done first, then measurements, then local corrections. And finally, we 
have inferred from this procedure that patterns with no dependencies, or us- 
ing only Pauli measurements, may only implement unitaries in the Clifford 
group. 

An obvious question is whether one can extend these ideas to other mea- 
surement based models, perhaps based on different families of entanglement 
operators, more general measurements and other types of local corrections. 
This is a matter which we wish to explore further. For now, it is already clear 
that both the notation and the calculus can be extended to the teleportation 
model which is based on 2-qubit measurements. This actually shows that 
teleportation models are embeddable in the one-way model in a very strong 
sense. We will return to this particular question elsewhere. 
We also feel that the methods explored here can be stretched further and 
made to be relevant to the study of error propagation and error correcting, 
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but this demands using mixed states, and interpreting patterns as cp-maps. 
Finally, there is also a clear reading of dependencies as classical communica- 
tions, while local corrections can be thought of as local quantum operations in 
a multipartite scenario. Along this reading, standardisation pushes non-local 
operations to the beginning of a distributed computation, and it seems the 
measurement calculus could prove useful in the area of quantum protocols. 
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9 Appendix 

We prove here that standardisation has indeed the properties quoted in the 
body of the paper. First, we need a lemma: 

Lemma 7 (Termination) For all ty, there exists finitely many *J}' such 

that y qy. 

Suppose ^3 has command sequence A n . . . Ai, and define for A t = d{AA = 
i, and for Aj = X*, d(Aj) = n — j. Define further: 

d(V) = (£ Ee ^(£),Ece^(C')) 
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This measure decreases lexicographically under rewriting, in other words 
*P =>- implies cZ(^3) > rf(^3'), where < is the lexicographic ordering on N 2 . 
Let us inspect all cases. First when one applies EC, then the first coordinate 
strictly diminishes (the second does not always, because of the duplication 
involved if C — X); when MC, the second strictly diminishes and the first 
stays the same or diminishes; when EA, the first strictly diminishes (because 
we dropped the case when A is itself an E), and maybe the second; when 
AC, the second strictly diminishes, and the first stays the same or diminishes 
(when A = E). 

Therefore, all rewritings are finite, and since the system is finitely branching 
(there are no more than to possible single step rewrites on a given sequence 
of length to), we get the statement of the theorem. □ 

It is not to difficult to strengthen the result above, by showing that the 
longest possible rewriting of *P is quadratic in to, where to is the length of ^'s 
command sequence. 

Say *}3 is standard if for no ^3', ^3 ^3'. 

Proposition 8 (Standardisation) For all ^J3, there exists a unique stan- 
dard ^3', such that *P =^>* *J3' ; and *J}' satisfies the (EMC) condition. 

Since the rewriting system is terminating, confluence follows from local con- 
fluence (meaning whenever two rewritings can be applied, one can rewrite 
further both transforms to a same third expression). Then, uniqueness of the 
standard form is an easy consequence (actually, for terminating rewriting sys- 
tems, unicity of standard forms and confluence are equivalent). Looking for 
critical pairs, that is occurrences of three successive commands where two 
rules can be applied simultaneously, one finds that there are only two types: 
EijMkCk with %, j and k all distinct, and E^M^Ci with k and I distinct. In 
both cases local confluence is easily verified. 

Suppose now *J3' does not satisfy (EMC). Then, either there is a pattern EA 
with A not of type E, or there is a pattern AC with A not of type C. In 
the former case, E and A must operate on overlapping qubits, else one may 
apply a free commutation rule, and A may not be a C since in this case 
one may apply an EC rewrite. The only remaining case in when A is of 
type M, overlapping E's qubits, but this is what condition (Dl) forbids, and 
since (Dl) is preserved under rewriting, this contradicts the assumption. The 
latter case is even simpler. □ 

9.1 Discussion 

This is what we wanted, namely we have shown that under rewriting any 
pattern can be put in (EMC) form. We actually proved a bit more, namely 
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that the standard form obtained is unique. 

However, one has to be a bit careful about the significance of this additional 
piece of information. Note first that unicity is obtained because we dropped 
the CC free commutations, and all EE commutations, thus having a very 
rigid notion of command sequence. One cannot put them back as rewrit- 
ing rules, since they obviously ruin termination and uniqueness of standard 
forms. 

A reasonable thing to do, would be to take this set of equations as gener- 
ating an equivalence relation on command sequences, call it =, and hope to 
strengthen the results obtained so far, by proving that all reachable standard 
forms are equivalent. 

But this is too naive a strategy, since Ei 2 X 1 X 2 = E 12 X 2 X 1 , and: 

Z? ys vt vs rys yt yt T? 

— VS 7t rys yt T? 
= vV 1 ^-^2^-2-^12 

obtaining an expression which is not symmetric in 1 and 2. To conclude, one 
has to extend = to include the additional equivalence X\Z\ = Z[X[, which 
fortunately is sound since these two operators are equal up to a global phase. 
We conjecture that this enriched equivalence is preserved. 
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